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Competition between the Kondo effect and Ruderman-Kittel-Kasuya-Yosida interaction in the 
two-impurity Kondo problem can be phenomenologically described by the Rasul-Schlottmann spin 
model. We revisit this model from the quantum entanglement perspective by calculating both the 
inter-impurity entanglement and the local Kondo entanglement, the latter being the entanglement 
between a local magnetic impurity and its spatially nearby conduction electron. A groundstate 
phase diagram is derived and a discontinuous breakdown of the local Kondo entanglement is found 
at the singular point, associated concomitantly with a jump in the inter-impurity entanglement. 

An entanglement monogamy holds in the whole phase diagram. Our results identify the important 
role of the frustrated cross-coupling and demonstrate the local characteristic of the quantum phase 
transition in the two-impurity Kondo problem. The implications of these results for Kondo lattices 
and quantum information processing are also brieffy discussed. 
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I. INTRODUCTION 

The past decade has witnessed the growing power of 
quantum entanglements developed from the quantum in¬ 
formation in understanding the novel quantum states 
and quantum phase transitions in solids [ll-Q • Condensed 
matter systems involve spin-^ objects as natural qubits 
and various spin exchange interactions as sources of quan¬ 
tum correlations. The present paper will explore from the 
quantum entanglement perspective the variable Kondo 
effect in a two-impurity Kondo model (TIKM), an in¬ 
triguing problem involving two typical kinds of spin ex¬ 
changes. 

It is known that Kondo systems, consisting of both 
itinerant electrons and local magnetic moments, nat¬ 
urally involves the single-ion Kondo effect and the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. 
The competition between them plays a key role in cor¬ 
related systems ranging from diluted mangetic alloys to 
heavy fermion compounds [bl, 0- The issue has been in¬ 
tensively investigated within the TIKM[8| where, in ad¬ 
dition to the antiferromagnetic (AFM) Kondo coupling 
between a magnetic impurity (or local moment) and its 
spatially nearby conduction electron, the two magnetic 
impurities are also coupled due to the RKKY interac¬ 
tion. There are two stable fixed points: the strong Kondo 
coupling limit with each impurity spins being completely 
quenched by the Kondo effect, and the strong RKKY 
interaction limit with a different antiferromagnetic spin 
singlet d, 0. The later situation has usually an AFM 
ordered groundstate in the concentrated Kondo lattice 
case 0- The TIKM has been realized in nanoscale de¬ 
vices where the observed Kondo signature varies with 
tunable RKKY interaction [Til IHl- In addition to the 
conventional spin-density wave quantum phase transi¬ 
tion, the variation of Kondo effect and its competition 
with magnetic order in the Kondo lattice systems may 


result in a local type of quantum phase transition sig¬ 
naled by a breakdown of Kondo effect. 0 

Theoretically, the primary focus is on the intermediate 
regime where the single-ion Kondo and RKKY energy 
scales, represented by Tk and Jr respectively, are com¬ 
parable to each other so that the p^sical properties of 
the two stable fixed points crossover p[l^. Of particular 
interesting is the case when the crossover is sharpened 
leading to a phase transition. 

Early numerical renormalization group [oj and confor¬ 
mal field theory studies [llj on the particle-hole symmet¬ 
ric TIKM revealed an unstable interacting fixed point 
at a finite ratio Jr/Tr ~ 2.2 where non-Fermi liquid 
behaviors such as the divergent staggered susceptibil¬ 
ity and specific heat coefficient were observed. It was 
soon pointed out by Rasul and Schlottmann [Tfij that 
these intriguing features can be understood phenomeno¬ 
logically by an spin-only effective model involving two 
impurity spins Sa, Sr, and two conduction electron 
spins Sc (a), Sc(6), as shown in Figdl As a result of 
many-body process, the interaction terms in the corre¬ 
sponding Hamiltonian Eq. m emerge from the low en¬ 
ergy regime of the original TIKM ji E 0. The Tk 
describes the splitting between the Kondo singlet and 
the spin triplet states, while the cross-coupling K rep¬ 
resents the interaction-induced frustration. It is known 
that the critical point may be replaced by a crossover 
in the absence of particle-hole symmetry |17l-[^ . while 
with this symmetry a large degeneracy at the critical 
point is observed [l7|, l20|, l23| . Compatible with this ob¬ 
servation, the Rasul-Schlottmann (RS) spin model in¬ 
deed exhibits an enhanced degeneracy at a special point 
P : (K/Tk = 1, JrITk = 2) corresponding to the critical 
point [l^. 

On the other hand, Kondo effect or Kondo screening 
is conceptually associated with the notion of ” Kondo en¬ 
tanglement” i[3,0- The Kondo models with a sin¬ 
gle or two magnetic impurities have been investigated 
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from the quantum entanglement perspective [23431) . For 
the single impurity Kondo problem, Kondo screening im¬ 
plies formation of a Kondo singlet groundstatejUj. This 
is an entangled state consisting of the two parts, one 
is the local magnetic impurity, another is the rest of 
the whole system. The entanglement between the two 
parts defines the so-called single-impurity Kondo entan¬ 
glement (SIKE), a measure of Kondo entanglement typ¬ 
ically quantified by the von Neumann entropy [33|. The 
SIKE in a gapless bulk follows the well-known scaling 
law of the thermodynamic entropy at large distance over 
the coherent length of the Kondo cloud [25|, [13 • Eor 
the TIKM, several different impurity-related entangle¬ 
ments are considered. Similar to the SIKE, one could 
consider the two-impurity Kondo entanglement (TIKE), 
i.e., the entanglement between the two impurities and 
the rest of the system. The TIKE can be also mea¬ 
sured by the von Neumann entropy [Hi [33|- Another 
useful quantity is the inter-impurity entanglement (HE), 
i.e., the entanglement between the two local magnetic 
impurities [m, Hoj, [s^l- Such HE can be quantified by 
concurrence or negativity [Hi lHI- If been shown that 
the HE is non-zero when the RKKY interaction is at least 
several times larger than the Kondo energy scale [Hi [sj- 
However, this feature alone does not sufficiently guaran¬ 
tee a true phase transition [ItI [38| nor necessarily imply 
a full suppression of the Kondo effect. 

One should notice that the SIKE and TIKE quan¬ 
tified by the von Neumann entropy in Kondo systems 
with two- or more impurities usually mix the contribu¬ 
tions from conduction electrons and other impurities [s^. 
hence these quantities are not as distinct as in the single 
impurity Kondo problem. Therefore, an alternative mea¬ 
sure of the Kondo entanglement, capable of characteriz¬ 
ing the variation of Kondo screening across the quantum 
critical point in generic multi-impurity systems, is highly 
desirable [io). 

In this paper, we investigate a local Kondo entangle¬ 
ment (LKE), namely, the entanglement between a mag¬ 
netic impurity and its near-by conduction electron only. 
The definition of the LKE in generic Kondo systems is 
described in Appendix A. Generally, suppression of this 
quantity should imply the complete destruction of the 
Kondo effect, though its connection with the impurity 
quantum phase transition as in the generic TIKM re¬ 
mains to be clarified. As a concrete example, the HE 
and LKE in the RS model are evaluated on equal-footing, 
both quantified by the concurrence or negativity. An en¬ 
tanglement phase diagram is then obtained, exhibiting 
the crossover from the Kondo singlet phase (with non¬ 
zero LKE) to the inter-impurity AEM phase (with non¬ 
zero HE). A critical point P emerges along the strong 
frustration line K = Tk where both LKE and HE show 
sudden changes and the groundstate wavefunction shows 
a discontinuity. In the following sections, we shall present 
the results of the exact solutions of the RS model, the 
entanglement phase diagram, as well as an entanglement 
monogamy. We also briefly discuss implications of these 


results for generic Kondo lattices and quantum informa¬ 
tion processing. 

II. MODEL AND SOLUTIONS 

The RS spin model, presumably taking into account 
the effective many-body process, is a fixed-point Hamil¬ 
tonian of the TIKM. It is described by[li| 

Hrs = Tk[Sa • Sc{a) + Sb • Sc{b)] + JrSa • Sb 

^ K[SA-Ub)PSB-Ua)]. (I) 

Where, Sa and Sb denote the spin operators of the lo¬ 
cal moments at sites A, H, respectively. The electron 
spin density at the spatial site r = a, 6 is denoted by 
^c{r) = with Ca{r) being the 

annihilation operator of conduction electrons and a the 
Pauli matrices. In the present two-impurity Kondo prob¬ 
lem, only the spin degrees of freedom are relevant while 
the charge degrees of freedom are frozen. Thus Sc{a) and 
Sc{b) are treated as two electron spins localized at the 
sites r = a and b respectively. The interaction parame¬ 
ters include T^, the single-ion Kondo temperature, Jr^ 
the inter site RKKY interaction energy, and A, the cross 
Kondo coupling between the local moments and the elec¬ 
tron spins, as shown in Eig JT|T^ . All these interaction 
terms are spin-5'1/(2) invariant. The direct Kondo cou¬ 
pling takes place between the impurity and conduction 
electron in the assigned nearest neighboring sites (A, a) 
or {B^b). The cross Kondo coupling term {K) is induced 
effectively in the low energy limit of the original TIKM 
and K Tr when the even/odd channel symmetry is 
imposed. Eor our purpose we will assume that Jr and K 
are tunable independently with respect to Tr. Because 
Eq. (I) is invariant under a combination of the impurity 
permutation Sa ^ Sr and the exchange Tr ^ A, the 
regime with strong frustration corresponds to K = Tr. 
Hence we only need to consider 0 < KfTR < I. In the 
following we set Tr = 1 without losing generality. 

The eigenstates of Hrs can be classified into six cat¬ 
alogs: two singlets, three triplets, and one quintet, ac¬ 
cording to the decomposition of tensor representations of 
the SU (2) group: 2(g)2(g)2(g)2 = 2x103x305 
as listed in Appendix B. These eigenstates are con¬ 
structed based on a complete set of the conventional basis 
15^(A), s^{a); S^{B), s^{b)) and labeled by the total spin 
5, its z-component Sz, and the parity (with respect to 
permutations of the two local moments) [I 3 . Because we 
consider Tr^ Jr^ and K are all AEM, the groundstate is 
among the mixed states of the two singlets of even parity. 
They are denoted by (each up to a normalization factor) 

1^(0+)) = IM;n) + in;^l) (2) 

+ KP^{\n;n) + \U;U)) 
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level spacing (energy gap) A vanishes at a special point 
P, corresponding to K = Tk = so that the singlet 

1^(0/+)) is degenerate with |^(o+)) at this point. Mean¬ 
while, the odd triplet is degenerate with 1^(0+)) 

for P = I^Jr > 2. Hence the model indeed shows a 
strong frustration along the line K = 1 and exhibits an 
enlarged symmetry at Pin. The precise wavefunctions 
across the P point can be determined by taking either 
limits {K = 1, Jr = 2 — e) and {K = 1, Jr = 2 ^ e), 
e ^ 0+. It readily reveals a discontinuity in |T('o+)) as 
shown in Appendix B. 


FIG. 1: (color online) The cartoon picture of the RS model: 
the local magnetic moments and the conduction electron spins 
are denoted by big blue circled and small green circled arrows, 
respectively. Tk, Jr, and K represent the single-ion Kondo 
energy scale, the RKKY interaction, and the interaction in¬ 
duced cross Kondo coupling. Various entanglements among 
these spins are defined in the text. 


III. ENTANGLEMENT PHASE DIAGRAM 

Now, we start from the conventional Kondo impurity 
entanglement, i.e., the SIKE defined as the entanglement 
between a local moment, say S{A), and the rest of the 
system, denoted by A. It is measured by the van Norman 
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with P(o) = —JrI^ — K/2 — 1/2 and A = 
v/(l + K- Jfl)2 + 3(K-1)2. 

Another relevant low energy states are from the odd 
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entropy Srike = -Tr:^A){pimp{J^)^riPimpi^)}, where 
pimp{A) is the reduced density matrix pimp{A) = Tr^p. 
Here, p = |Tg')(Tg'| is the density matrix for the ground- 
state ( in our case = 1^(0+))) of the whole system. 
It is straightforwardly seen that Esike = 1 due to the 
5'P(2)-spin invariance, indicating a maximal entangle¬ 
ment between a local moment and the reminder of the 
whole system. 

Next, we consider the TIKE, the entanglement of 
the two local moments with the conduction electrons. 
Eollowing Ref. [29| , this entanglement is determined 
by the reduced density matrix of the two impurities, 
pimp{AB) = Tr(^c)Pi with indicating trace over the 

Hilbert subspace spanned by the conduction electrons. It 
can be quantified by the von Neumann entropy (2^ 


Etike = -Ps logps-{l-Ps) log 


1 -Ts 


( 8 ) 


(Ut;W-|^t;W) 


Here, Ps = \—fAB is the fidelity of the spin singlet within 
the reduced two impurity state, and /ab = (^o+|>5(A) • 
5'(P)|To+) is the spin-spin correlation function on the 
groundstate. In our case, 

1 l^K -Jr 1 

~ 4‘ 
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^{jR-K-iy^3{K-iy 
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is then evaluated on the groundstate |To+) 
as shown in Eigl2l We find that Etike is not only a 
smoothly varying function of Jab as already shown in 
Ref. [2^, but also a smooth function of K and Jr with¬ 
out detectable feature across the point P. 

(6) Now we turn to the HE. It can be measured by 
the concurrence or negativity , Cue- Its evaluation is 
also related to the reduced two-impurity density matrix 
pimp{AB). The concurrence can be expressed by [2^ 


(7) 


Cue = max{—2 /ab — 1/2,0}. 


( 10 ) 


It is apparent that the groundstate of the RS Hamilto¬ 
nian Eq.([T]) is the singlet |T('o+)). However, the repulsive 


The result is plotted in Eigl3}a). Eor fixed K < 1, Cue 
increases continuously with Jr. Eor K = 1, Cue shows a 
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FIG. 2: (color online) The TIKE as a function of Jr and K. 


sudden increase from zero to unity when Jr goes across 
Jr = 2. 

Together with the observed discontinuity in the 
groundstate wavefunction, the jump in HE evidences a 
transition at P [i^ . But its relation with the suppression 
of Kondo effect remains uncertain. We now turn to an al¬ 
ternative definition of the Kondo entanglement, i.e., the 
LKE between a local moment, say S{A), and the conduc¬ 
tion electron at its nearest neighbor site, Sc{a). The LKE 
differs to the conventional impurity entanglement as it 
involves only a spatially neighboring pair formed by a lo¬ 
cal moment and a conduction electron. This local Kondo 
pair is defined in the original real space point-contact 
Kondo interaction JrSa * Sc{a) (with Jr being the orig¬ 
inal Kondo coupling) as shown in Appendix A. Similar 
to HE, the LKE can be evaluated by the concurrence or 
negativity, via the corresponding reduced density matrix 
pLR{Aa) = T'^(^Ad)P^ with indicating the trace in 

the Hilbert space except the subspace spanned by Sa and 
Sc{a). Thus we have 


Clre = max{—2/^a — 1/2, 0}, (11) 


where fAa = (^o+|5'(A) • Sc(a)|4^o+) is the correlation 
function of the local Kondo singlet state. 


1 JrP2K -A 1 

4 ^{Jn -K-lf + 3{K - if " 4‘ 


( 12 ) 


The result of Crre is plotted in FiglSj^b). Interestingly, 
Clre develops a maximum for 0 < iF < 1, < 2 and 

decreases monotonically for Jr > 2. But along the line 
K = 1 it shows a sudden suppression Jr > 2 . 

An entanglement phase diagram in terms of K and 
Jr is then drawn in Fig 01 where three different phases 
divided by the lines Jr = K P 1 and Jr = A — 2K are 
indicated: the HE phase {Cue > ^fihRE = 0), the LKE 
phase {Clre > 0,C//£; = 0), and the co-existence phase 
{Cue > 0, Clre > 0). The HE and LKE phases contact 
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(b) 

EIG. 3: (color online) The concurrence of the HE (a) and the 
LKE (b) as functions of Jr and K. The single-ion Kondo 
energy scale is taken as unit Tr = 1. 


only at the point P: by increasing Jr across P along the 
strong frustration line iF = 1, /ab has a sudden drop 
from 1/4 to —1/2, /ag has a jump from —1/2 to 0. Or, 
Cue has a jump from 0 to 1 while Clre has a sudden 
drop from 1/2 to 0. 

Therefore, together with the discontinuity of the wave- 
function 14^(0+)), the sudden changes along the line K = 
1 in the HE and LKE do evidence a phase transition ac¬ 
companied by a breakdown of Kondo effect. Of course, 
a true second order phase transition usually involves a 
continuous variation of the order parameter before its 
suppression. So the discontinuity exhibited in the HE or 
LKE ( as an order parameter here) is seemingly due to 
the simplicity of the present model involving only two 
conduction electrons. 

We notice that the inter-impurity spin-spin correlation 
function, closely related to the HE discussed here, was ac¬ 
tually calculated in the early numerical renormalization 
group study @ at low yet finite temperatures, where the 
calculated quantity did not exhibit a sudden discontinu¬ 
ity but a rather sharp change at the critical point. Inter- 
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FIG. 4: (color online) Entanglement phase diagram: there 
are three distinct phases divided by the lines Clke = 0 and 
Cue = 0. These two lines intersect at the singular point 
P : (K = 1, = 2), where the coexistence regime diminishes 

leading to a Kondo breakdown transition to the AFM singlet 
phase. 


estingly, a more recent calculation based on the natural 
orbitals renormlization group method does find a suddent 
jump of this quantity at zero tempertureQ- Therefore, 
the discontinuity of the HE and LKE may be not limited 
to the present model, but a generic feature of the impu¬ 
rity state involving finite degrees of freedom at quantum 
critical points. This discontinuity could be smeared when 
the impurity degrees of freedom become infinite as in the 
Kondo lattice case. 


IV. THE FIDELITY OF THE KONDO AND 
INTER-IMPURITY AFM SINGLETS 

In order to clarify whether the phase with non-zero 
Cue or Clke > 0 corresponds to the AFM inter¬ 
impurity or Kondo singlets, we calculate the normalized 
wavefunction overlaps and re¬ 

spectively. We consider the pure state of the Kondo 
screening phase (denoted by as the groundstate 

at the fixed point K = Jr = {). Similarly, we de¬ 
note \^aem) the pure inter-impurity AFM state at 
the fixed point A = 0, Ji? ^ oo. Here, = 

|5'(0)Aa) ® \'^{0)Bb), I'^AFm) = |4'(0)ylB) ® |5'(0)ab), with 
l^(O)Aa) = ^l)-| W), l^(O)Bb) = -^,i\ ^l)-l ^t)), 

i^(o)ab) = M)-i W), i^(o)ab) = -^{\ t;)-i ;t)). 

We find that in the regime with Cue = 0 or Clke = 0 
the respective wavefunction overlap vanishes. Therefore, 
the obtained entanglement phase diagram Fig H] reflects 
the overall evolution from the Kondo singlet to the inter¬ 
impurity AFM singlet. 

Specifically, increasing K will enhance the frustration, 
so that the fidelities of and aem) on the true 

groundstate |^g) change with varying Jr. The fidelities 
can be calculated as the normalized wavefunction over- 
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(b) 

FIG. 5: (color online) (a) The fidelity of the Kondo singlet 
on the groundstate |Tg). (b) The fidelity of the inter¬ 
impurity AFM singlet |Tafm) on the groundstate |Tg). 


laps and The three-dimensional 

plots of the wavefunctions overlaps are shown in Fig|5]^a) 
and Fig. ETb), respectively. FigEfa) and FigEl^b) show 
the Ji^-dependence of these quantities for fixed values 
of K. We find that (^gI^afm) increases rapidly with 
Jr and saturates to 1 after entering the AFM phase, 
Jr > 2. Notice that (^gI^ks') always approaches to 
1/2 for Jr oc, indicating that the Kondo singlet fi¬ 
delity Ps ^ 1/2. As shown in Appendix A, this implies 
the non-correlated Kondo state or a full suppression of 
Kondo effect. 


V. ENTANGLEMENT SUM RULE 

By definition, the concurrences of LKE and HE intro¬ 
duced in the main text. Cue and Clke, are always non¬ 
zero ( but < 1) when the respective correlation function 
/ab < —1/4 or fAa < —1/4. Based on the exact solu- 
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(a) 




FIG. 6: (color online) The groundstate fidelity of the Kondo 
singlet ITks') (a) or of the inter-impurity AFM singlet (b) as a 
function of the RKKY interaction for various fixed frustration 

K. 


tions, these functions can be reexpressed as 

fAB = -^cosO- = ^cos6>- ^ sin6> - ^.(13) 

Where, 9 = tan-^ ^ = \jIZk-i\ - Therefore, 

in the whole parameter regime Jr > 0, K < 1, or0< 
0 < f, we have 

Cue = max{Acos6>, 0}; (14) 

A a/3 

Clke = max{--cos6>+— sin6>,0}. (15) 

Thus Cue and Crke are both non-zero when K 1 < 
Jr < 4 — 2K. In this co-existence regime where ^ < 0 < 
we have 

Cue + Clke = sin(6> + —). (16) 

6 

This is a sum rule constraining the variations of Cu e and 
Clke- ^ < Cue + Clke < 1 in this regime. 


As an inference of the sum rule, we note that a Werner 
state with Ps ^ (1 + 3/v/2)/4 0.78 has the nonlo¬ 

cal correlation characteristic, i.e., violation of the Bell 
inequality [m |45|, |46|. This value corresponds to C ^ 
0.56. Therefore, we have an entanglement monogamy: 
Clke and Cue cannot simultaneously maximize, nor 
even simultaneously fall into this region. In other words, 
only one of the entanglements could be maximized or 
violate the Bell inequality. 

Finally, one can also define the entanglement between 
the spin Sa and Sc(6), i.e., the cross Kondo entanglement, 
quantified by Ccke- Owing to the discrete symmetry 
mentioned previously, this quantity can be drived similar 
to CLKE , obtaining Ccke = max{ — ^ cos 0 — ^ sin 0}. 
So CcKE is non-zero only when Jr < AK — 2. This region 
is a subregion of the LKE phase, below the straight line 
connecting {K = 0.5, Ji? = 0) and the P-point(not shown 
in FiglU). Hence including Coke does not violate the 
previous sum rule. Moreover, another sum rule Clke + 
CcKE = cos^ holds in this subregion. Therefore the 
previous inference as well as all the conclusions in the 
main text remain valid in the whole phase diagram. 


VI. IMPLICATIONS AND DISCUSSIONS 

As we emphasized, the quantum entanglement related 
to the single-ion Kondo effect is complicated when two 
or more local magnetic impurities are introduced. We 
have considered the LKE in the case of two magnetic 
impurities where a quantum phase transition from the 
Kondo to the AEM singlets takes place. The exact re¬ 
sults obtained from a simplified yet effective TIKM, i.e., 
the RS spin model, show a phase diagram involving sev¬ 
eral regimes corresponding to non zero LKE and HE, and 
their co-existence. With increasing frustrated cross cou¬ 
pling the co-existence regime shrinks and vanishes at the 
critical point P where the groundstate has a discon¬ 
tinuity accompanied by jumps in the HE and LKE. The 
discontinuity of such quantities at the critical point may 
be a generic feature caused by the non-extensive impurity 
term in the free energy of the TIKM [dll li^. 

Due to the entanglement sum rule, the LKE and 
HE cannot be simultaneously maximized even in the 
co-existence regime. Because an entangled state with 
non-zero concurrence while still keeping the Bell-CHSH 
inequality [m |46| can be used for quantum information 
processing including quantum teleportation [13, [56|, the 
singlet groundstate in the co-existence phase shares this 
feature either by the LKE or HE, and thus should be 
an interesting candidate state for quantum information 
processing. 

It is interesting to understand implications of our 
present study for more generic Kondo lattice models. In 
these systems, like the heavy fermion metals, the mag¬ 
netic quantum phase transitions may be influenced by 
the variation of Kondo effect. A local quantum phase 
transition is possible in more generic Kondo lattice phase 
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diagram where the criticality is associated with a critical 
breakdown of the single-ion Kondo effect jiol-lsH . Near 
the critical point, the Hall constant shows a discontin¬ 
uous jump due to the reconstruction of the Fermi sur¬ 
face across the critical point [i^, [5^. Experimental ev¬ 
idence for this scenario comes from several prototypes 
of heavy fermion metals including YbRh2Si2[52| and 
CeNiAsO[53| where the observed Hall constant exhibits 
a sudden change accompanying the magnetic phase tran¬ 
sition. 

Although the above scenario could be naturally un¬ 
derstood based on the Kondo entanglement picture, a 
lattice model Hamiltonian with exact solutions showing 
the Kondo entanglement breakdown transition is lack¬ 
ing. Thus our studied model can serve as a toy model 
to understand its basic physics from the quantum entan¬ 
glement perspective. On the one hand, one expects that 
with the evolution from two impurities to an regular lo¬ 
cal moment lattice, the Kondo singlet and inter-impurity 
AFM singlet states evolve into the paramagnetic heavy 
fermion and AFM ordered phases, respectively. On the 
other hand, in addition to the RKKY interaction, the 
crossing Kondo coupling {K) emerges as a many-body 
frustration effect and plays a role in controlling the tran¬ 
sition. Generally, in the regime with relatively small 
K, the Kondo and inter-impurity singlets can co-exist in 
the intermediate regime of Jr. The co-existence regime 
diminishes with increasing K and a direct Kondo sin¬ 
glet breakdown transition takes place when the Kondo 
coupling is maximally frustrated (at K = Tr). In the 
Kondo lattice phase diagram this condition should cor¬ 
respond to the regime with strong geometric frustrations 
and quantum fluctuations but no spin liquid phase sets 
in[5l|. Finally, extending the present study to the Kondo 
lattice cases with frustrated cross couplings or various 
anisotropic interactions is not an easy task but highly de¬ 
sirable. Nevertheless, our present result already provides 
a concrete example of Kondo breakdown quantum phase 
transition from the quantum entanglement perspective. 
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Appendix A: The LKE in the A-impurity Kondo 
system 

The Kondo system consisting of a metallic electron 
host, with energy dispersion e(/c), and A-number of quan¬ 
tum mechanical magnetic moments (of spin-1/2) local¬ 


ized at sites 2,2 = 1,2,.,A, is described by the Hamil¬ 

tonian 

Hnkm = Y / e{k)c\{k)ca{k) 

N 

+JKY.Si-Uri), (Al) 

i=l 

where the Kondo coupling is either antiferromagnetic 
{Jk > 0) and locals in the sense that it takes place 
between a local moment Si and a spatially nearby con¬ 
duction electron at site fi via the point-like interaction 
JxSi • Sc(ri), with the electron spin operator Sc{ri) = 
^ct{ri)daa'Ca'{ri). 

The local Kondo entanglement (LKE) refers to the 
quantum entanglement of a local Kondo pair consisting 
of an impurity fixed spin Si and its spatially nearby con¬ 
duction electron. As we focus on the spin degrees of 
freedom relevant in the Kondo screening, such a pair of 
spins is dubbed as a local Kondo state. It constitutes a 
subsystem, with a local Hilbert space fti spanned by Si 
and Sc{ri). Let O be the total Hilbert space, and the 
Hilbert subspace complementary to Qr. QiUQi = Q. As¬ 
sume I^g) be the groundstate of the A-impurity Kondo 
model Hr km defined in Eq. (jAl[) . p = |Tg')(4^g| the cor¬ 
responding density matrix of the whole system, the LKE 
is defined as the entanglement of the reduced density ma¬ 
trix of this subsystem obtained by taking trace over all 
other degrees of freedom except the local Hilbert space: 
pQ. = Tr^ p. Obviously, has dimension 4, so pQ. can 
be expressed by the matrices / 4 x 4 , r 0 a, with r being 
the Pauli matrices of the local spin. Si = Owing 

to the fact that the Hamiltonian is real, SU{2) invariant, 
and the groundstate is always a spin-singlet, one has the 
general form [3, 

pn, = 1/4x4 + J[rr ® at + rf 0 af + rf ® af], (A2) 
In terms of the Bell basis, the maximal entangled states 

= ;^(l ^;>±l W) and |$(±)) = ^(1 ^t)±l W)), 

the reduced density matrix can be expressed by 

pQi = (A3) 

with Ps and pt being the probabilities of spin singlet and 
spin-triplet, respectively. Thus, pQ. is a mixture of spin- 
singlet and spin-triplet. Ps and pt are related to the spin- 
spin correlation function fs = •s(ri)|TG') viaps = 

\ — fs^Pt = \ Therefore, the pure local Kondo 

singlet corresponds to Ps = 1 or while the 

pure local Kondo triplet corresponds to pt = -^ or fs = 
Notice that the case of | or p^ = Pt = 1/2 

corresponds to the non-correlated Kondo pair with equal 
mixture of spin-singlet and spin triplet. 
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In our definition, the LKE is measured by the con¬ 
currence of pq .. Because . can be reexpressed 
as pn, = ^^|^(-))(^(-)| + ^i^/ 4 x 4 , it is also a 
Werner state with the spin singlet fidelity Ps = 
The the concurrence of this state is 

given bvj^ 


Appendix B: Eigenstates and eigen energies of the 
Rasul-Schlottmann model 

The sixteen eigenstates (each up to a normalization 
factor) and the corresponding eigenvalues are solved as 
following: 


Clke = max{2ps — 1,0}. (A4) 

Such entanglement can be also measured by the 
negativity [361 of the reduced density matrix, JVp^.. It is 
straightforwardly seen that the negativity of the Werner 
state is equal to the concurrence [^. i.e., N'p^. = Clke- 

Similarly, the entanglement of a pair of two local mag¬ 
netic moments, say Si and Sj^ can be defined following 
the previous approach. Parallel to the LKE, such the 
inter-impurity entanglement (HE) measured by the con¬ 
currence {Cue) is closely related to the spin-spin cor¬ 
relation function • Sj\^G)- This quantity varies 

monotonically with the RKKY interaction jR{ij). The 
indirect RKKY interaction, i.e., the coupling between the 
local moments which is not explicitly present in the bare 
Hamiltonian, is generated in the second order perturb- 
tion in Jk and dependend on the density of states at 
the Eermi energy and the spatial separation | Ri — Rj \ of 
the two local moments. Because JR{ij) oscillates non- 
universally, the IIE’s of different pairs (i,i) are usually 
complicated. Remarkably, there are two special situa¬ 
tions where the HE competes with the LKE: (i) the two 
impurity case with N = 2; (ii) The sites i and j are 
the nearest-neighboring sites in the Kondo lattice with 
N = L {L = the total number of sites). Therefore, 
the competition between the LKE and HE manifests the 
competition between the Kondo singlet and the inter¬ 
impurity AEM singlet in the TIKM, or manifests the 
competition between the paramagnetic heavy Eermion 
state and the AEM ordered state in the Kondo lattices. 

A phenomenological understanding of above competi¬ 
tion invokes two energy scales, the single-ion Kondo tem¬ 
perature Tr ^ D^/pFJK and the RKKY 

interaction Jr^ with D the band width and pp the den¬ 
sity of states at the Eermi energy. In the renormaliza¬ 
tion group treatment one approaches the low energy limit 
by gradually integrating out the degrees of freedom of 
conduction electrons. So that at T = 0, the latter can 
be introduced by adding a direct RKKY term jRSi • Sj 
into the original Hamiltonian, while the single-ion Kondo 
screening effect is described by a term TpSi • Sc{ri). 
Phenomenologically, Tr is the excitation energy of the 
Kondo triplet above the Kondo singlet groundstate. The 
effective cross Kondo coupling (denoted by K in the main 
text) between Sj and Sc{fi) can be also induced by pure 
quantum mechanical many-body processes. Proper inter¬ 
impurity distance and the direct RKKY interaction guar¬ 
antee the required particle-hole symmetry [I^- These in¬ 
terpretations provide a basis for the RS spin model as a 
minimal fixed point Hamiltonian of the TIKM. 


1^(0+)) = I W;n) + l^t;^;) + ci(|M;W + IM;M)) 
-(ci + i)(|^t;^l) + IM;W) (Bi) 


E, 


( 0 +) 


Jr 
' 4 


K 1 

y ~ 2 


^{1 + K- JrY + 3(K - 1)2 


(B2) 


IV+)) = IM;^t) + in;^;) + c2(|M;W + in;n)) 

-(c2 + i)(|M;tl) + IM;^t)) (B3) 


E(o'+) — 


Jr_E 

4 2 

1 , y(i+if-Ji?)2+3(i^-i)2 


(B4) 


|v|/(i+,+i)) = I n) - I ^t) + I tt; W - I 

0)) = (B5) 


-B(i+) = 


Jr K 1 


(B6) 


+c 4 (|M;n)-|M;W) 

14/(1-,_i)) = C3(in;^;)-|^;;M)) 


(B7) 


^ _ Jr ViK - 1)^ + 4 

yi-) --R -^- 


(B8) 


IV-.+d) = C5(|M;n)-|^t;W) 
IV-.O)) = 


(B9) 
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+ 

1 

II 

1 


-1) 

2 

2 _|_ 72 


(BIO) 

1^(2, 

2)) = 







1^(2, 

d) = 

l^t;n) + l 


+ 1 


+ 


1^(2, 

0)) = 

1 W;n) + I 


+ 1 

^t;n) 

+ 




ffi hi) T 




(Bll) 

^(2,- 

-d) = 

ftl; hi) + 


+ 1 


+ 


^(2.- 

-2)) = 

IU;U) 







cq = ^ ^ ^ Along the line K = 1 oi 

Jr = 2K^ some coefficients are divergent, but the cor¬ 
rect forms can be obtained by taking the limits from ei¬ 
ther sides. In particular, in the vicinity of the singular 
point P : (iC = 1, Ji? = 2), the wavefunctions are de¬ 
termined unambiguously by taking the limits at P ap¬ 
proached from either sides along the line K = 1. For 
instance, when P = 1, Ji? = 2 + e, e ^ 0+, we have the 
normalized state 

1^(0+))+ = 

+l^t;M) + IM;W]- (B13) 
While when P = 1, = 2 — e, e ^ 0+, we have 


E(2+) = 


Jr K I 


(B12) 


In above, ci 


v/(l+K-Jfi)2+3(K-l)2-JH+2 


-^{i+K-jRy+KK-iy-jR+2 
Jr-2K 


Jr-2K 
C3 = 


C2 = 

JR+^J{K-1Y+Jl 
1-K ’ 


C4 = 


1-K+V(K-1)^ + Jg 


Jr 


C5 = 


JR-^y{K-ly+J% 


1-K 


and 


l^'(o+))- = ^[2|n;^;) + 2iM;^t)-in;W) 

(bi4) 

Therefore, +(T('o+) |T('o+))- = 0. This result demon¬ 
strates a discontinuity of the groundstate wavefunction 
across the singular point. 


[1] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 
416, 608 (2002). 

[2] T.J. Osborne and M.A. Nielson, Phys. Rev. A 66, 032110 

( 2002 ). 

[3] G. Vidal, J.I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. 
Lett. 90, 227902 (2003). 

[4] L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. 
Mod. Phys. 80, 517 (2008). 

[5] R. Horodecki, P. Horodecki, M. Horodecki, and K. 
Horodecki, Rev. Mod. Phys. 81, 865 (2009). 

[6] A.G. Hewson, The Kondo Problem to Heavy Fermions^ 
Gambridge University Press, Gambridge, England, 1993. 

[7] P. Goleman, in Handbook of Magnetism and Advaneed 
Magnetie Materials^ V.l, 95 (Wiley, 2007). 

[8] G. Jayaprakash, H.R. Krishnamurthy, and J.W. Wilkins, 
Phys. Rev. Lett. 47, 737(1981). 

[9] B.A. Jones and G.M. Varma, Phys. Rev. Lett. 58, 843 
(1987); B.A. Jones, G.M. Varma, and J.W. Wilkins, 
Phys. Rev. Lett. 61, 125 (1988). 

[10] S. Doniach, Physica B+G 91, 231 (1977). 

[11] L.I. Glazman and R.G. Ashoori, Science 304, 524 (2004). 

[12] N.J. Graig et al.. Science 304, 565(2004). 

[13] P. Gegenwart, Q. Si, and E. Steglich, Nature Phys. 4, 
186 (2008). 

[14] B.A. Jones and G.M. Varma, Phys. Rev. B 40, 324 
(1989); B.A. Jones, B.G. Kotliar, and A.J. Millis, Phys. 
Rev. B 39, 3415 (1989). 

[15] 1. Affleck and A.W.W. Ludwig, Phys. Rev. Lett. 68, 1046 
(1992) 

[16] J.W. Rasul and P. Schlottmann, Phys. Rev. Lett. 62, 
1701 (1989); P. Schlottmann and J.W. Rasul, Physica B 
163, 544 (1990). 

[17] 1. Affleck, A.W.W. Ludwig, and B.A. Jones, Phys. Rev. 
B 52, 9528 (1995). 


[18] O. Sakai, Y. Shimizu, and T. Kasuya, Solid State Gom- 
mun. 75, 81 (1990). 

[19] R.M. Eye, Phys. Rev. Lett. 72, 916 (1994). 

[20] J. Gan, Phys. Rev. Lett. 74, 2583 (1995). 

[21] J.B. Silva et al., Phys. Rev. Lett. 76, 275 (1996). 

[22] G. Zarand, G.H. Ghung, P. Simon, and M. Vojta, Phys. 
Rev. Lett. 97, 166802 (2006). 

[23] There are two types of particle-hole symmetry in the 
TIKM with lattice inversion symmetry, depending on 
the inter-impurity distance being even or odd lattice 
spacing, respectively. Usually, only the first type with 
even inter-impurity distance is the required symme¬ 
try which may gua rantee a critical point with ground- 
state degeneracy [3, [ij. As an effective spin-only model, 
Eq.([T]) does not show such distinction for even/odd inter¬ 
impurity distance. Instead, the required groundstate de¬ 
generacy can arise by tuning the cross Kondo coupling, 
see the following discussions. 

[24] A.T. Gosta, Jr. and S. Bose, Phys. Rev. Lett. 87, 277901 

( 2001 ). 

[25] K. Le Hur, P. Doucet-Beaupre, and W. Hofstetter, Phys. 
Rev. Lett. 99, 126801 (2007). 

[26] K. Le Hur, Ann. Phys. 323, 2008 (2008). 

[27] 1. Affleck, N. Laflorencie, and E.S. Sorensen, J. Phys. A 
42, 504009 (2009). 

[28] A.T. Gosta, Jr., S. Bose, and Y. Omar, Phys. Rev. Lett. 
96, 230501 (2006). 

[29] S.Y. Gho and R.H. Mckanzie, Phys. Rev. A 73, 012109 
(2006). 

[30] E. Eriksson and H. Jonhannesson, Phys. Rev. B 84, 
041107(R) (2011). 

[31] H. Saleur, P. Schmitteckert, and R. Vasseur, Phys. Rev. 
B 88, 085413 (2013). 

[32] K. Yosida, Phys. Rev. 147, 223 (1966). 


















10 


[33] C.H. Bennetti, Phys. Rev. A 54, 3824 (1996). 

[34] A. Bayat, S. Bose, P. Sodano, and H. Johannesson, Phys. 
Rev. Lett. 109, 066403 (2012). 

[35] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[36] G. Vidal and R.F. Werner, Phys. Rev. A 65, 032314 

( 2002 ). 

[37] M.-F. Yang, Phys. Rev. A 71, 030302(R)(2005). 

[38] The discontinuity in the first derivative of negativ¬ 
ity usually comes from the requirement of nonnegative 
concurrence [ 3 ^ . Its relation with a true quantum phase 
transition should be further corroborated by checking 
the nonanalyticity of groundstate energy or groundstate 
wavefunction. 

[39] In some circumstances such as in psudo-gaped or dissi¬ 
pative hosts the SIKE measured by the von Neumann 
entropy can show a peak feature at the quantum critical 
points [ 2 I, [ 2 ^. 

[40] In a spin chain model with two magnetic impurites sev¬ 
eral impurity-related entanglement quantities measured 
by the negativity display detectable features in their 
derivatives at the quantum phase transition [3. These 
quantities do not exhibit apparant variation of the Kondo 
effect across the transition. More recently, the Schmidt 
gap has been proposed as an order parameter measuring 
the impurity quantum phase transitionsfiH]. 

[41] A. Bayat, H. Johannesson, S. Bose, and P. Sodano, Nat. 
Commun. 5, 3784 (2014). 


[42] In the conformal field theory and non-Abelian bosoniza- 

tion approaches the enlarged hidden SO (7) symmetry at 
the non-trivial critical point is identified [id. . 

[43] L.-A. Wu, M.S. Sarandy, and D.A. Lidar, Phys. Rev. 
Lett. 93, 250404 (2004). 

[44] R.-Q. He, J. Dai, and Z.-Y. Lu, arXiv:1501.01834, 2015. 

[45] J.F. Clauser, M.A. Horne, A. Shimony, and R.A. Holt, 
Phys. Rev. Lett. 23, 880 (1969). 

[46] R. Horodecki, P. Horodecki, and M. Horodecki, Phys. 
Lett. A 200 , 340 (1995). 

[47] S. Popescu, Phys. Rev. Lett. 72, 797 (1994). 

[48] J. Lee and M.S. Kim, Phys. Rev. Lett. 84, 4236 (2000). 

[49] Q. Si, S. Rabello, K. Ingersent, and J.L. Smith, Nature 
413, 804 (2001). 

[50] P. Coleman, C. Pepin, Q. Si and R. Ramazashvili, J. 
Phys. C 13, R723 (2001). 

[51] Q. Si, Physica B: Conden. Mat. 378, 23 (2006). 

[52] S. Paschen et ah. Nature 432, 881 (2004). 

[53] Y.K. Luo et ah, Nature Materialsl3, 777 (2014). 

[54] M. Nielson and 1. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[55] R.F. Werner, Phys. Rev. A 40, 4277 (1989). 

[56] S. Lee, D.P. Chi, S.D. Oh, and J. Kim, Phys. Rev. A 68, 
062304 (2003). 


